Nucleon-deuteron capture with chiral potentials 
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Present day chiral nucleon-nucleon potentials up to next-to-next-to- 
next-to leading order and three nucleon forces at next-to-next-to leading 
order are used to analyze nucleon-deuteron radiative capture at deuteron 
laboratory energies below ~ 100 MeV. The differential cross section and 
the deuteron analyzing powers Ay (d) and Ayy are presented and compared 
to data. The theoretical predictions are obtained in the momentum-space 
Faddeev approach using the nuclear electromagnetic current operator with 
exchange currents introduced via the Siegert theorem. The chiral forces 
provide the same quality of data description as a combination of the two- 
nucleon AV18 and the three-nucleon Urbana IX interactions. However, the 
different parametrizations of the chiral potentials lead to broad bands of 
predictions. 

PACS numbers: 21.45. -Hv, 25.10.-f-s, 25.40.Lw 

1. Introduction 

Various modern (semi)phenomenological nucleon-nucleon (NN) poten- 
tials 121 El adjusted to the nucleon-nucleon data allow for a relatively 
good description of few-nucleon systems. Moreover, the existing models 
of a three nucleon force (3NF) when combined with those potentials led 
generally to an improvement in the description of three and four nucleon 
bound states HHSHnHIl and reactions in the 3N-continuum |HliniEl- De- 
spite indisputable successes such forces lack reliable theoretical justification 
in the underlying theory of strong interactions - the quantum chromody- 
namics (QCD). Recently, a new family of nuclear potentials has been de- 
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rived |1H I12| ITH l I14j which has a direct connection to QCD. Incorporating 
the spontaneously and exphcitly broken chiral symmetry, which plays a fun- 
damental role in QCD, it is possible to construct nuclear forces in the frame- 
work of chiral perturbation theory. This approach has been originally pro- 
posed by Weinberg further developed by Ordonez et al. [TH], Kaiser 
et al. J7j, Entem et al. and by Epelbaum et al. (See also recent 
review articles |191 1121 [TT] and references therein). Up to now, NN and 3N 
forces have been worked out and applied in few-nucleon systems up to next- 
to-next-to-next-to-leading order (N^LO) and next-to-next-to-leading order 
(N^LO), respectively. The leading 4N- forces have also been worked out re- 
cently 120]. The chiral forces have already been successfully used to describe 
the elastic nucleon-deuteron scattering 1^21^ and the deuteron breakup pro- 
cesses \21\ I22j . We remind the reader that applications of chiral nuclear 
forces are limited to the low-energy region. 

Chiral effective field theory has also been used to study processes with 
external pions and/or photons, see [221 IMl 123 US] for some recent applica- 
tions. While the neutron-proton radiative capture has already been consid- 
ered in this framework T7| , no such analysis is available for electromagnetic 
processes with three nucleons. In this study we would like, at least partially, 
to fill this gap by presenting for the first time the chiral predictions for the 
n + d ^ 'J +^ H and p + d ^ 7 +^ He processes. 

We use the chiral NN potentials jl4j at next-to leading order (NLO), 
next-to-next-to leading order (N^LO) and at next-to-next-to-next-to lead- 
ing order (N^LO) of chiral expansion. In the N^LO where the 3NF con- 
tributes for the first time, we also include such 3NF's flTI consistent with 
the two-body interaction at this order. A complete study of electromagnetic 
processes in the chiral effective field theory would require the usage of the 
nuclear electromagnetic current operator derived consistently with the ap- 
plied nuclear forces. Such a current operator is, however, not yet available, 
see Ref. j28| for some early work along these lines. For that reason, we 
decided to adopt in the present study our standard approach ^\ , in which 
the electromagnetic current operator is taken in the form of a one-body 
current supplemented with some many-body contributions introduced via 
the Siegert theorem. 

In Section 2 we shortly describe our theoretical formalism. Results and 
a comparison to data are presented in Section 3. We summarize in Section 
4. 

2. Theoretical formalism 

The present study is based on the formalism described in detail in PHIIHUI 
I31j . The calculations are done for the two-body photodisintegration of the 
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3N bound state and the nuclear matrix element for the radiative A^d-capture 
N™'^ is found by applying the time reversal to the photodisintegration am- 
plitude A^^'^. The latt er one is a matrix element of the electromagnetic 
current operator between the initial 3N bound state | ^'h) and the final 

scattering state | ^'^j) 

^ (^y I J, I ^^.) • (1) 

In the Faddeev scheme A'^^'^ can be written in the form j3()j 

iV^^^=(Vi|(l + P)jJ^b), (2) 

where | is a Faddeev component of the state | ^jvi) -P is a permu- 
tation operator defined as a sum of cyclical and anti-cyclical permutations 
of three particles 

P = P12P23 + P13P23, (3) 
where Pij interchanges the i-th and j-th nucleons. 

Starting from the Schrodinger equation for the scattering state | ^'yvJ) 
the transition amplitude can be written as 

n;^'' = {^i I (1 + p) I I + {<Pi\p\ u), (4) 

where | (pi) is a product of the internal deuteron state and a momentum 
eigenstate of the free nucleon-deuteron motion. The auxiliary state | U) 
fulfills the Faddeev-like equation 

1 if) = (^tGo + ^(1 + P)vl'^Go{tGo + 1)) (1 + P)j^, \ ^b) 

+ (^tGoP + ^(1 + P)v}'^Go{tGo + I U), (5) 

where v}^^ is a part of the 3NF which is symmetrical under the exchange 
of nucleons 2 and 3, Gq is the free three-nucleon propagator, and ti is the 
two-body t-operator acting in the 2-3 subspace. 

Equation ^ is solved in the momentum-space basis 

\p,q,aj) =\p,q,ils)j,{X,^)I{jI)J;{t^)T), (6) 

where p and q are the magnitudes of the two Jacobi momenta. The quantum 
numbers I, s, and j are the two-body subsystem orbital angular momentum, 
spin and the total angular momentum, respectively. The orbital angular 
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momentum A of the spectator nucleon together with its spin couples to 
the total spectator angular momentum I. The angular momenta j and I 
couple finally to the total angular momentum of the 3N system J. Similar 
coupling {t^)T occurs for the isospin quantum numbers. 

Starting from the transition matrix element A'^^'^ the observables follow 
in a standard way pIT] . 

3. Results 

The chiral potentials depend on low energy constants (LEG) and the 
cut-off parameters A and A. The LEC's are determined for given A and 
A from a fit to NN and 3N data. The A cut-off appears in the spectral 
function regularization and ensures that the high-momentum compo- 
nents of the two-pion exchange are explicitly excluded from the potential. 
The cut-off A comes from a further regularization of the potential as used 
in the Lippmann-Schwinger equation, which cuts off the (meaningless) con- 
tributions of high-momentum states in the dynamical equation. We per- 
formed numerical calculations for a few combinations of (A, A) values which 
covers the possible ranges of the regularization parameters: (450(400 at 
NLO), 500), (600(550 at NLO),500(600 at N^LO)), (550,600), (450(400 at 
NLO),700), (600(550 at NLO),700) (in units of [MeV/c]). 

The theoretical predictions at a given chiral order provide a band com- 
prising all results. The bands presented in the following correspond to the 
NN potential at NLO, N^LO, NN-F3NF at N^LO, and finally to the NN 
potential at N'^LO. 

Let us start the presentation of our results with the cm. differential cross 
section for the neutron-deuteron (Figs.nnj and proton-deuteron (Figs. HIS)) 
capture. In Fig. ^ we display predictions at incoming neutron lab. energy 
En=9 MeV and 10.8 MeV. The light (cyan) band corresponds to NLO pre- 
dictions while the dark (red) one refers to the full (i.e. NN-I-3NF) N^LO 
results. While the NLO band is relatively broad, the N^LO band is much 
narrower. This shows the fast convergence of the chiral expansion for this 
observable. The NN-I-3NF N^LO chiral predictions are in agreement with 
the AV18-|-Urbana IX combination of the standard NN and 3NF models 
obtained using the same current operator. In Fig.[21we show a scale of 3NF 
effects at N^LO by comparison predictions which corresponds to the full 
Hamiltonian (dark band) and other ones which corresponds to the inter- 
action restricted to the NN force only (light band). Results obtained with 
NN force at N^LO are similar to those at NLO. Due to the width of the 
NN band it is difficult to conclude if the inclusion of the 3NF improves or 
not the description of the data, however the 3NF effects are clear. There is 
an improvement at smaller Q^n angles but at larger angles the data favour 
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the pure NN predictions. The comparison between N^LO and N^LO based 
on NN interactions only is shown in Fig. |31 The inclusion of the additional 
terms appearing in N^LO (without corresponding 3NFs) shifts the predic- 
tions above those of N^LO and brings them above the data. It would be 
interesting to see if the inclusion of the 3NFs in N^LO will decrease the 
cross section bringing theory back to the data, similarly as it is the case at 
N^LO. 

In Figs. 14161 proton-deuteron capture cross sections are shown at two 
incoming deuteron lab. energies, E£;=29 and 95 MeV. The nice convergence 
of chiral expansion is seen at both energies (see Fig. In Fig. |21 at smaller 
deuteron energy Ed=29 MeV the inclusion of the 3NF in N^LO shifts the 
results in the wrong direction, deteriorating the description of data. Pass- 
ing to (incomplete) N^LO improves the agreement with data. At the much 
higher deuteron energy Erf=95 MeV the 3NF acts in the opposite direction, 
undoubtedly improving the data description. At this energy there is no 
difference between predictions at N^LO and N'^LO based on the NN inter- 
actions only. However, the band of predictions at N^LO is much broader 
than the corresponding one at N^LO. 

Now let us discuss selected spin observables. We would like to show 
results for the deuteron vector analyzing power Ay{d) fFigs. 171^ and the 
deuteron tensor analyzing power Ayy (Figs. 110112(1 at four deuteron energies 
17.5 MeV, 29 MeV, 45 MeV and 95 MeV, for which comparison with data 
is possible. Here, the convergence of the chiral expansion is much slower 
(see Figs. [7| and llUp and the bands have more or less the same width. At 
the highest deuteron energy Erf=95 MeV for Ay{d) the N^LO band is even 
wider than that one in NLO. The effects of the 3NF at N^LO on Ay{d) are 
not substantial and are seen mainly at the highest deuteron energy. The 
observed discrepancy with the data might at least partially be caused by 
the simplified model of the electromagnetic operator used in the present 
calculations. (The results based on the AV18-|-Urbana IX potential and a 
more realistic model of the current operator is shown in 31 where a better 
agreement with data was found.) The effects of the new terms occurring in 
the NN potential in N^LO are small, as seen in Fig. |^ At the three lower 
energies they slightly shift predictions in the data direction. At the highest 
energy the picture for Ay{d) is similar when including the 3NF or the new 
N^LO terms. However, in the first case the band of predictions is much 
broader. 

Similarly to the Ay{d) case, the effects of moving from NLO to N^LO, 
as well as the 3NF effects for Ayy (Figs. I1UII11() are not significant. The 
discrepancy between AV18-|-Urbana IX and the chiral predictions for the 
higher energies of 45 MeV and 95 MeV is particularly interesting. At these 
energies the step to the next order in the chiral expansion seems to be 



6 



pdcap'chiral060614 printed on February 9, 2008 



necessary (see Fig. I12|) . Moreover, the additional terms in the N^LO NN 
chiral potential shift predictions farther away from the data. This points to 
the need of developing the corresponding 3N forces, which is under way. 

4. Summary 

The first predictions for radiative nucleon-deuteron capture based on 
chiral forces are presented. The differential cross section and the deuteron 
vector and tensor analyzing powers have been calculated using the chiral 
NN potential at NLO, N^LO and N^LO. Also the results at N^LO with the 
chiral 3NF of that order are shown. The nuclear current operator contains, 
besides the single nucleon current, some many-body contributions as intro- 
duced via the Siegert theorem. The chiral potentials give similar predictions 
as the AV18+Urbana IX combination of the standard nuclear forces. We 
found, however, a significant theoretical uncertainty for all considered ob- 
servables. The bands of predictions are in some cases too broad, especially 
in the N^LO NN and N^LO NN+3NF cases, to make a decisive conclusion 
about the quality of the Nd-capture data description. This broadness can 
be partially caused by the inconsistent current operator used in this study. 
In future studies the nuclear current constructed in the same framework as 
chiral forces should be used. This first study indicates that the inclusion 
of the consistent current operator and going to N3L0 in the chiral expan- 
sion might be necessary to obtain a more precise description of radiative 
processes within the chiral effective field theory framework. However, we 
have shown that already now this framework can be used for analyzing low 
energy electromagnetic processes. 
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Fig. 1. The cm. neutron-deuteron capture cross section at: 9.0 MeV (a) and 10.8 
MeV (b) neutron lab energies. The hght (cyan) and dark (red) shaded bands are 
NN NLO and NN+3NF N^LO predictions, respectively. The solid line represents 
prediction obtained with the AV18 and Urbana IX forces. Data are from |32 |. 



8 



pdcap'chiral060614 printed on February 9, 2008 




Fig. 2. The cm. neutron-deuteron capture cross section at: 9.0 MeV (a) and 10.8 
MeV (b) neutron lab energies. The hght (cyan) and dark (red) shaded bands are 
NN and NN+3NF force predictions at N^LO. The hue and data as in Fig.Q 
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Fig. 3. The cm. neutron-deuteron capture cross section at: 9.0 MeV (a) and 10.8 
MeV (b) neutron lab energies. The hght (cyan) and dark (red) shaded bands are 
NN force predictions at N^LO and N'^LO, respectively The solid line and data as 
in Fig.n 
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Fig. 4. The cm. proton-deuteron capture cross section at: 29 MeV (a) and 95.0 
MeV (b) deuteron lab energies. The bands and the hne as in Fig. ^ Data at 29 
MeV are from and at 95 MeV from 
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Fig. 5. The cm. proton-deuteron capture cross sections at: 29 MeV (a) and 95 
MeV (b) deuteron lab. energies. The bands and the line as in Fig. |21 Data as in 
Fig.H 
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Fig. 6. The cm. proton-deuteron capture cross sections at: 29 MeV (a) and 95 
MeV (b) deuteron lab. energies. The bands and the line as in Fig. O Data as in 
Fig.H The wider N^LO band covers the N^LO at Ed=95 MeV. 



pdcap 'chiraW60614 printed on February 9, 2008 



13 




Fig. 7. The deuteron vector analyzing power Ay{d) at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as 
in Fig.n Data at 17.5 MeV are from jSSlEHj, at 29 and 45 MeV from |37|. at 95 
MeV from 




Fig. 8. The deuteron vector analyzing power Ay{d) at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as in 
Fig. Data as in Fig. □ 
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Fig. 9. The deuteron vector analyzing power Ay{d) at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as in 
Fig. Data as in Fig. □ 



16 



pdcap'chiral060614 printed on February 9, 2008 




Fig. 10. The deuteron tensor analyzing power Ayy at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as in 
Fig.n Data at 17.5 MeV from at 29 MeV from |3H] (squares) and jST] 

(circles), at 45 MeV from 1371 (circles) and (squares) and at 95 MeV from |34| . 




Fig. 11. The deuteron tensor analyzing power Ayy at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as in 
Fig. Data as in Fig. UHl 




Fig. 12. The deuteron tensor analyzing power Ayy at: 17.5 MeV (a), 29 MeV (b), 
45 MeV (c) and 95 MeV (d) deuteron lab. energies. The bands and the line as in 
Fig. Data as in Fig. UHl 
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